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1 Introduction 

In Q and the subsequent paper Connes gave a rather intrinsic construction of a Hilbert 
space intimately associated with the zeros of the Riemann zeta function on the critical line. 
But the zeros having multiplicities higher than a certain level (which is a parameter in Connes's 
construction), have (if they at all exist) their contributions limited to that level, and not to the 
extent given by their natural multiplicities. Thus subsists the problem of a natural definition 
of a so-called "Hilbert-Polya space", with orthonormal basis indexed by the zeros p of ^ and 
integers k varying from to rup — 1 where nip is the multiplicity of p. We do not solve that 
problem here but we do propose a rather natural construction of Hilbert space vectors X^j^, 
C{p) = 0, k < nip, which in the limit when the parameter A goes to become perpendicular 
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(when they correspond to distinct zeros. The vectors corresponding to a multiple root p are 
independent but need to be orthogonalized.) As in Connes's constructions these vectors live in 
a quotient space. Controlling the limit A — > to obtain a so-called Hilbert-Polya space probably 
involves considerations from mathematical scattering theory (we have previously studied in |^ , 
1^] some connections with the problems of L— functions.) 



The context in which our construction takes place is that of the Nyman-Beurling formulation 
of the Riemann Hypothesis as an approximation problem |14], [^]. Let K = L^(]0, cx3[, dt) 



(over the complex numbers), let x be the indicator function of the interval ]0, 1], and let p be 
the function "fractional part" (the letter p is also used to refer to a zero of the Riemann zeta 
function, hopefully no confusion will arise.) Let < A < 1 and let B\ be the sub-vector space 
of K consisting of the finite linear combinations of the functions 1 1— > p{j)-, for A < ^ < 1. 



Theorem 1.1 (Nyman [14], Beurling Q) The Riemann Hypothesis holds if and only if 



XG U 

0<A<1 



Actually we are following Q here in using a slight variant of the original Nyman-Beurling 
formulation. It is a disappointing fact that this theorem can be proven without leading to any 
new information whatsoever on the zeros lying on the critical line (basically what is at works is 
the factorization of functions belonging to the Hardy space of a half-plane ||l^.) The following 
is thus rather remarkable: 



Theorem 1.2 (Baez-Duarte, Balazard, Landreau and Saias Q) Let us write D{X) for 
the Hilbert-space distance infjgg^ \\x — f\\- We have 
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If the Riemann Hypothesis fails this result is true but trivial as the left-hand side then takes the 
value +00. So we will assume that the Riemann Hypothesis holds. The sum on the right-hand 

side is over all non-trivial zeros p of the zeta function, counted only once independently of their 
multiplicities nip. We prove the following: 



Theorem 1.3 We have: 




So the zeros are counted according to the square of their multiplicities. To prove this lower 
bound we will construct remarkable Hilbert space vectors X'^^ , ^(p) = 0, fc < nip and use them 
to control D{X). The following "toy-model" gives us reasons to expect that the lower bound in 
fact gives the exact order of decrease of D(A): 



Theorem 1.4 Let Q{z) = Yl^^i^ — a • z)'^" be a polynomial of degree q > 1 will all its roots a 
on the unit circle (the root a having multiplicity vria). Let P{z) be an arbitrary polynomial. Let 



E{N,P):= inf / \P{z) - Q{z)A{z)\ 
We have as N goes to infinity: 

lim NE{N,P) =^ml \P{a)\ 



2^ 
27r 



2 The prediction error for a singular MA(q) 



As motivation for our result we first consider a simpler approximation problem, in the context 
of the Hardy space of the unit disc rather than the Hardy space of a half-plane. Let Q{z) = 
^0.(1 — a • z)"^" be a polynomial of degree q> 1 will all its roots a on the unit circle (the root 
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a having multiplicity nia so that q = Xlo'^a-) define: 

E{N):= inf / \1 - Q(z)A{z)\' - 
deg{A)<N J SI 27r 

The measure ^ is the rotation invariant probability measure on the circle S^, with z = exp(i0). 
The minimum is taken over all complex polynomials A{z) with degree at most A^. We are 
guaranteed that limjv^oo ^(-^) = as Q{z) is an outer factor (ll^). More precisely: 



Theorem 2.1 As N goes to infinity we have: 

lim NE{N) =^ml 

a 



Note 2.2 In case Q{z) has a root in the open unit disc then E{N) is bounded below by a 
positive constant. In case Q{z) has all its roots outside the open unit disc, then the result 
above holds but only the roots on the unit circle contribute. Finally if all its roots are outside 
the closed unit disc then the decrease is exponential: E{N) = 0{c^), with c < 1. 



The theorem, although not stated explicitely there, is easily extracted from the work of Grenan- 
der and Rosenblatt They state an 0(;^) result, in a much wider set-up than the one 

considered here (which is limited to simple-minded q-th order moving averages.) Unfortunately 
the bound is now believed not to be systematically true under their hypotheses (as is 

explained in |^]; I thank Professor W. Van Assche for pointing out this fact to me.) Never- 
theless their technique of proof goes through smoothly in the case at hand and yields the exact 
asymptotic result as stated above. We only sketch briefly the idea, as nothing beyond the tools 
used in |]l^] is needed. 

We point out in passing that it is of course possible to express E{N) explicitely in terms of the 
Toeplitz determinants for the measure dfi = |Q(exp(i6'))|2|^. But already for an MA{2) this 
gives rise to unwieldy computations. . . . Rather: let be the vector space of polynomials of 
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degrees at most N + q, let Vn be the subspace of polynomials divisible by Q{z), and let Wn be 
its g-dimensional orthogonal complement. Then E(N) is the squared norm of the orthogonal 
projection of the constant function 1 to Wn- A spanning set in Wn is readily identified: to 
each root a one associates Y^q, Y^-^, . . . , ^^^a-i defined as 

Y^o:=l + az + ... + a''+iz''+i 

yj^i := z + 2az^ + . . . + (AT + g)a^+9-i^Af+</ 

and similarly for k = 2, . . . , rria — 1- We can then express E{N) using a Gram formula in terms 
of (the inverse) of the positive matrix (of fixed size q >^ q but depending on N) built with the 
scalar products of the y's. It turns out that in the limit when N goes to infinity and after 
the rescaling i— > X^^^ := N^^^'^^'^Y^^ the Gram matrix decomposes into Cauchy blocks 
(l/(i + j + l))o<jj<ma of size ma, one for each root a. It is known from Cauchy that the 
top-left element of the inverse matrix is m^. This is how ^^J"" arises, after keeping track 
of the scalar products (1,X^^). Instead of the constant polynomial 1 we could have looked 
at the approximation rate to an arbitrary polynomial P{z). The proof just sketched applies 



identically and gives the Theorem 1.4 from the Introduction 



3 Invariant analysis and a construction of Baez-Duarte 

The Mellin transform f{t) ^ /(s) = /^^q f{t)t'^^^dt establishes the Plancherel isometry between 
K = L2(]0,oo[,(it) and L'^{s = ^ + ir, with inverse F{s) ^ L=i/2+iT^(^)*~'i^^- ^et a{s) 
be a measurable function of s (as a rule when using the letter s we implicitely assume Re(s) = ^. 
We will use letters w and z for general complex numbers.) If a(s) is essentially bounded then 
F{s) I— > a{s)F{s) defines a bounded operator on K which commutes with the unitary group 
De '■ f{t) ^ "Tl-^^i)' bounded operators commuting with the Dg {0 < 9 < oo) are 

obtained in such a manner. More generally all closed invariant operators are associated to a 
measurable multiplier a{s) (finite almost everywhere, but not necessarily essentially bounded). 
For the details of this technical statement, see 0]. 
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For example the Hardy averaging operator M : f (t) t-^ j ^-^ f{u)du corresponds to the spectral 
multiplier The operator 1 — M corresponds to the spectral multiplier a-^d is thus 

unitary. Another (see Q) remarkable invariant operator is the (even) "Gamma" operator 
r+ = !F-^I. Here / is the inversion f{t) i— s- jf{j) and .7-+ is the additive Fourier transform as 
applied to even functions (the cosine transform). The multiplier associated to r_|_ is the (Tate) 
function 

/ \ i-s r(s/2) C(l - s) _5 r'^s N N /"°° sin(27ru) 

7+ s = 7r2 ^ / ; = =2^ 'tt ' cos — r s = 1 - s / ^ ^ du 

r((l-s)/2) C(s) 2 7o 

A further invariant operator is the operator U introduced by Baez-Duarte Q in connection 
with the Nyman-Beurling formulation of the Riemann Hypothesis: its spectral multiplier is 
so [/ = (M - = J^+I{M - 1). 

From the results recalled above on invariant operators, we see that invariant orthogonal projec- 
tors correspond to indicator functions of measurable sets on the critical line. So a function f[t) 
is such that its multiplicative translates D0{f ) (0 < 9 < oo) span K if and only if F{s) = f{s) 
is almost everywhere non-vanishing (Wiener's L^-Tauberian Theorem.) In that case the phase 
function 

is almost everywhere defined and of modulus 1. It thus corresponds to an invariant unitary 
operator, also denoted Uf. 



Let us introduce the anti-unitary "time-reversal" operator J acting on K as g ^id)- The 
operator Uf commutes with the contractions-dilations, is unitary, and sends / to J(/). We call 
this the Bdez-Duarte construction as it appears in |jl| (up to some non-essential differences) 
in relation with the Nyman-Beurling problem (the phase function arises in other contexts, 
especially in scattering theory.) 

To relate this with the operator U = {M — l)jr_^/, one needs the formula 

as) 



Jo ^^i^^"'^* 
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which is fundamental in the Nyman-Beurhng context. This formula shows that U is the phase 
operator associated with p{j)- 

Generally speaking, the operators Uf are related to the Hardy spaces M2 = L^{]0,l],dt) and 
H^"'" = oo[, dt) (we will also use the notation for the Mellin transform of L'^{]0, l],dt).) 

Indeed the time-reversal J is an isometry (anti-unitary) between and H^"'". Let us assume 
that the function / belongs to H^. The operator Uj has the same effect as J on f, but contrarily 
to J is an invariant operator. This puts the space B\{f) (of finite linear combinations of 
contractions Dg{f) for X < < 1) isometrically in a new light as a subspace of L^([A, oo[, dt). 
The marvelous thing is that in this new incarnation it appears to be sometimes possible to find 
vectors orthogonal to B\{f) and thus to get some control on B\{f) as A decreases (as in the 
Grenander-Rosenblatt method.) 



4 The vectors Yg\ 



To get started on this we first replace the function — with an element of H^. This is 
elementary: 



Proposition 4.1 (pi, [|T0|) The function Z{s) = belongs to H^. Its inverse Mellin 



transform A{t) is given by the formula 



^(t) = [^]log(t) + log([l]!) + [i] 



One has (forO< t < 1) A{t) = i log(i) + 0(1). 



The Baez-Duarte construction will then associate to A{t) the operator V with spectral multiplier 

V{s)- 



s V c(i-^) 



l-sj C{s) 
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so that 

V = {l- M f ■ U 

This last representation will prove useful as it allows to use the formulae related to U from 
1^] and IQ]. Let (0 < A < 1) be the sub- vector space of H-^ of linear combinations of the 
contractions Dai A) for X < 9 < 1. The function = - — i- is the Mellin transform of 

Xlit) '■= (1 + log(t))x(0- The quantity D{X) considered by Baez-Duarte, Balazard, Landreau 
and Saias is thus the Hilbert space distance between Xi{t) To bound it from below we 

will exhibit remarkable Hilbert space vectors X^j, indexed by the zeros of the Riemann zeta 
function and perpendicular to C\. We then compute the exact asymptotics of the orthogonal 
projection of xi to the vector spaces spanned by the X^^^, for a finite set of roots, exactly as in 
the Grenander-Rosenblatt method. 

To each complex number w and natural integer > we associate the funtion 'tpw,k{t) = 
(log( j))^ t^"^ x{t) on ]0, oo[. For Re(t(;) < 1 it is integrable, for Re(zi;) < ^ it is in K. Let Q\ be 
the orthogonal projector from K onto L'^{[X, oo[). The main point of this paper is the following: 

Theorem and Definition 4.2 For each < A < 1, each s on the critical line, and each 
integer k > the L'^-limit in K of V^'^Q\V{il)w,k) exists as w tends to s from the left half- 
plane: 

:= U.m. V-^QxV{^^,k) 

Rc{m)<| 

For each X <9 <1 the scalar products between Dq{A) and the vectors Y^j^ are: 

Note 4.3 The proof shows the existence of an analytic continuation in w accross the critical 
line, but we shall not make use of this fact. 

Clearly one has the following statement as an immediate consequence: 



9 



Corollary 4.4 Let < A < 1. The vector Y^f, is perpendicular to C\ if and only if C^\s) = 
for all j < k, if and only if s is a zero p of the zeta function and k < rup. 

Note 4.5 Our scalar products (/, 5) are complex linear in the first factor and conjugate- linear 
in the second factor. 

Note 4.6 The operator ^ when applied to a not necessarily analytic function on the critical 
line is defined to act as (where s = ^ + ir.) 

Proof The proof of existence will be given later. Here we check the statement involving the 
scalar product, assuming existence. The following holds for A < < 1 and Re(u;) < \: 

{V-^QxV{i^y,,k),De{A)) = {QxV{i^^,k),VDe{A)) 

= {QxV{i^.,,k),De-V{A)) 

= {Vi4^^,k),Qx-D0-J{A)) 

= {V{tl;^,k),De-V{A)) 

= {tlJyj,k , De{A)) 

= (^) i^'^'~'"i~'"xm,A) 

\dwj 7]o,i] 

Taking the limit when w ^ s gives 

= f--r-l [ t-^-'^A{{)dt 

\tdTj J]o,i] 
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Taking the complex conjugate: 

k 



De{A),Y^^,) = (.A^ 6^- l^^^^t--2+^^ A(t)dt 

[ A{t)dt 



k 

ds) 



^ \ k 

d \ „„ 1 



ds 

which completes the proof (assuming existence 



To prove the existence we will use in an essential manner the key Lemme 6 from Q. We 
have seen that = (1 — M)'^U where M is the Hardy averaging operator and U the Baez- 
Duarte operator. The spectral function U{s) extends to an analytic function U{w) in the strip 
< Re(w) < 1. We need pointwise expressions for V{ipw^k){t), t > (at first only Re(ty) < ^ 
is allowed here). Thanks to the general study of U given in [P, we know that for Re(u') < ^ 
the vector U{tpw^k) in K is given as the following limit in square mean: 

d sin(27rt/?;) 



1- f r^ /InnA- -w d sm(27rt/?;) , 
l.i.m. / (log(-))''t; ^ / ' dv 



Following with a slight change of notation, we now study for each complex number w with 
Re(?i;) < 1 (and each integer A: > 0) the pointwise limit as a function of t > for 5^0: 

/,N 1- f^n -w d sin(2TTt/v) , 

ip^,k(.t) := ]im / log - S \ / ' dv 

Jg V dv irt/v 



Theorem 4.7 ([§]) Let k = 0. For each t > and Iie{w) < 1 the pointwise limit defining 
^w,o{t) exists. It is holomorphic in w for each fixed t. When w is restricted to a compact set in 
Re(tt;) < 1, one has uniformly in w the bound 92^1, o(^) = 0{j) on [1, oo[. Uniformly with respect 
to w satisfying < Re(?x)) < 1 one has (pwfiit) = U {w) t~^ + 0(1) on < t < 1. 



Proof Everything is either stated explicitely in Lemme 6 and Lemme 4, or follows from 
their proofs. We will give more details for A; > 1 as this is not treated in • 
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Corollary 4.8 For each w in the critical strip < Re(u;) < 1 the Hardy operator M : f{t) — > 
i /o ■^(^) applied arbitrarily many times to ipw,o{t)- The functions M^{(pyj^o) (X G 

are ) 077, [1^ oo[, uniformly with respect to w when it is restricted to a compact subset 

of the open strip, and satisfy on t &]0,1] the estimate M^{(pu),Q)it) = (i^TTy) U{w)t~'^ + 0(1), 
uniformly with respect to w. 

Proof A simple recurrence. • 
We thus obtain: 

Corollary 4.9 The vectors Y^^ exist (for Re(s) = One has the estimates: 

m^)(t) = o(ii±M^))!) (,e[i,ooD 

V{Y,>^^om = V{s)t-' + 0{l) {X<t<l) 
^(^^o)(*) = (0 < i < A) 
uniformly with respect to s when its imaginary part is bounded. 

Theorem 4.10 Let k > 1. For each t > and Re(u;) < 1 the pointwise limit defining fw,k{t) 
exists. It is holomorphic in w for each fixed t. When w is restricted to a compact set in 
Re(it;) < 1, one has uniformly in w the bound ^w,k{t) = on [l,oo[. Uniformly for 

< Re(u;) < 1 one has (pw,ki^) = (f^H i""') + 0{1) onO<t<l. 

Proof The formula defining ipw,k{t) is equivalent to (after integration by parts and the change 
of variable u = 1/v): 

^^,k{t)= lim i- r(A: + ii;logH)(log(tx))^-^ isin(27rin) — 
A-»oo7riJi ^ 
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This proves the existence of ipw,k{t), its analytic character in w, and the uniform 0{j) bound 
on [1, oo[. The formula can be rewritten as: 

f d \ w ,„ 1 / N du 

Vw,k{t) =[-r] u'^-^smilTTtu) — 

\dw J irt Ji u 

When w is in the critical strip the integral u"^^^ sin{2TTt u) ^ is absolutely convergent and 
its value is t^"^ J^^ u""'^ sm{2iTu) = j^{2tt 1^-"" cos{^)r{w) from well-known integral 
formulae, so that: 

k 



I — w 2 T^t Jn u J 



d 

dw ) \1 — w 2 ' ^ ' T^i Jo ^ ' ) 

The first term is (^)'^ (^(^) i""") and the second term can be explicitely evaluated using the 
series expansion of sin(27rt u) with the final result 

which shows = (^)*' (C/(w) t""') + 0(1), on < t < 1, uniformly for < Re(ii;) < 1. • 



As was the case for /c = we then deduce that the Hardy operator can be applied arbitrarily 

A 

I, A: 



many times to i-Pw.k for < Re(?i;) < 1. The existence of the follows. 



Theorem 4.11 Let k > 0. The vectors Y^j, exist (for Re(s) = One has the estimates: 

viy.^km = o(il±M!))!) (t.[i,oc[) 

ViY^'km = {Vis)t-n+0{l) (A<t<l) 

ViYs^m = (0<t<A) 

the implied constants are independent of A and are uniform with respect to s when its imaginary 
part is bounded. 



Proof Clearly a corollary to 4.1C. 
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5 The vectors and completion of the proof 

Definition 5.1 Let < A < 1. To each zero p of the Riemann zeta function on the critical 
line, of multiplicity nip, and each integer < k < rrip we associate the Hilbert space vector 

/ \ — — ~k 

where Y^^ = 1 .i.m.^ — V ^ Q {'4^w,k^ > ^ the unitary operator (il<f — I , Qx is orthogonal 
projection to L'^{[X,oo[,dt), and il^w^kit) = (log( j))^ t""" x(*)- 



Note 5.2 Of course there is no reason except psychological to allow only zeros of the Riemann 
zeta function at this stage. 



Theorem 5.3 As A decreases to one has: 



lim = (pi^P2) 

A — >(J 

lim = i 

A^o ^ P^^' P^^' k + l + l 



Proof To establish this we first consider, for Re(si) = Re(s2) = \: 



If si 7^ S2 an integration by parts shows that it is O (log(-^)) ^ ^. On the other hand when si 

1 



S2 its exact value is -^j^-^j^^ (log(^))'''^^''^^"^. With this information the theorem follows directly 



from 4.11 as (for example) the leading divergent contribution as A —> to {y(Xsk)^ ^O^sS) 
is V{s)V(J)!l (if t~' t-^dt which gives (log(i)) The rescaling Y ^ X 

is chosen so that a finite limit for {X^^,X^^) is obtained. As the scalar products involving 
distinct zeros have a smaller divergency, the rescaling let them converge to 0. • 
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Theorem 5.4 Let Xii't) = (1 + log(t))x(i)- As A decreases to one has: 



1^ ^ ^.A 



lim yiog(-) (xi,X;^fc) = {k>l) 



Jim A/log(i)(xi,X,^.) ^ ^ 



Proof We have (1 - M)xi = X, and V = {l- Mf U so Vxi = (1 - M) C/x- From || we know 



that Ux is so yxi is the function 



1 f't sin(27rt)) 



1 r* 
t Jo 



dv. It is thus O(t^) as t ^ 0, and 



from 4.11 we then deduce that the scalar products (xii^A) admit finite hmits as A ^ 0. This 



settles the case k>l. For /c = 0, one uses the uniformity with respect to u) in 4.7 to get 



lim (xi , Yp o) = Mm (xi , iPwfi) 
which gives lim^_^p (1 + log(t)) f^'"^ dt = 



p-i 



We can now conclude the proof of our estimate. 



Theorem 5.5 We have: 



liminfD(A)yiog(-)> 



E 



Proof Let i? be a non-empty finite set of zeros. We showed that D{\) is the Hilbert space 
distance from xi to C\, and that the vectors X'^^ for < /c < rup are perpendicular to C\. So 
D(A) is bounded below by the norm of the orthogonal projection of xi to the finite-dimensional 
vector space Hr spanned by the vectors X'^j^, < k < nip, p ^ R. This is given by a well-known 
formula involving the inverse of the Gram matrix of the X\^s as well as the scalar products 



{xi,X^f^). From 5.3 the Gram matrix converges to diagonal blocks, one for each zero, given 



by Cauchy matrices of sizes nip x nip. From Cauchy we know that the top- left element of the 



inverse matrix is trip. Combining this with the scalar products evaluated in |5.4| we get that the 
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squared norm of the orthogonal projection of xi to Hr is asymptotically equivalent as A ^ 
to —T—TTT—- The proof is complete. • 



We can apply our strategy to a fully singular MA(q) on the unit circle. The relevant Baez- 
Duarte phase operator will then be (up to a non-important constant of modulus 1) the operator 
of multiplication by z^'^ and it is apparent that this leads to a proof equivalent to the one we gave 
in our previous discussion, inspired by |11]. In the case of the Nyman-Beurling approximation 



problem for the zeta funtion, we expect in the quotient of by C\ a "continuous spectrum" 



additionally to the "discrete spectrum" provided by the (projection to of the) X\ 's, C,{p) 



0, /c < rup. It is tempting to speculate that the rescaling will kill this continuous part as A — > 0, 
so that in the end only subsists a so-called "Hilbert-Polya" space. This would appear to require 
|5.5| to give the exact order of decrease of the quantity D{\) and the numerical explorations 
reported by Baez-Duarte, Balazard, Landreau and Saias in Q seem to support this. 
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